A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 
02.08.1; 08.09.3, 08.18.1. 



ASTRONOMY 

AND 

ASTROPHYSICS 

1.2.2008 



o 

Q 



> 

o. 

(N- 

Oh; 
6' 

cd ■ 



% 



Mixing processes during the evolution of red giants with 
moderate metal deficiencies : the role of molecular-weight 
barriersB 

Corinne Charbonnel^'^, JefFery A. Brown'^'^, George Wallerstein^'^ 

^ Laboratoire d'Astropliysique de Toulouse, CNRS UMR 5572, Toulouse, France 

^ Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 

^ Dept. of Astronomy, University of Washington, Box 351580, Seattle, WA 98195-1580, USA 

■* Program in Astronomy, Washington State University, Pullman, WA 99164-3113, USA 

^ Guest observer. Dominion Astrophysical Observatory, Herzberg Institute of Astrophysics, National Research Council, Canada 

Corinne.Charbonnel@obs-mip.fr, jbrown@delta.math.wsu.edu, wall@astro.washington.edu 

February 1, 2008 



Abstract. We have assembled accurate abundance data 
for Li, C, and N as well as the ^^C/^^C ratio for five field 
giants with [Fe/H] ~ —0.6 including Arcturus and two 
stars in the globular cluster 47 Tuc. Using their very pre- 
cise Mboi values obtained from HIPPARCOS parallaxes, 
we can place them into an evolutionary sequence. The se- 
quence shows that the '^^C/^^C ratios drops from ~ 20 
to near 7 between Mboi = +1 and -1-0.5, while Li disap- 
pears. At the same time the ^^C/^'*N ratio diminishes by 
0.2 to 0.4 dex. The two stars in 47 Tuc with Mboi near 
-2.0 show even lower ^^C/^'*N ratios by 0.4 dex indicating 
further mixing as they evolved to the top of the red giant 
branch. 

These observations confirm the existence of an extra- 
mixing process that becomes efficient on the red giant 
branch only when the low-mass stars reach the so-called 
luminosity function bump. We use the values of the carbon 
isotopic ratio observed in our sample to get constraints on 
the /j,-barriers that may shield the central regions of a star 
from extra-mixing. We show that the same value of the 
critical gradient of molecular weight leads to ^^C/^'^C ra- 
tios observed at different metallicities. This "observational 
critical /x-gradient" is in very good agreement with the one 
which is expected to stabilize meridional circulation. This 
result provides strong clues on the nature of the extra- 
mixing which occurs on the RGB, and indicates that it is 
related to rotation. 
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1. Introduction 

For some time now it has been clear that evolved low mass 
stars exhibit chemical anomalies which are not predicted 
by standard stellar evolution theory. By "standard" we 
refer to the hypotheses that stellar convective regions are 
instantly mixed and that no transport of chemicals oc- 
curs in the radiative regions. In addition, standard mod- 
els are those without rotation at any depth. With these 
assumptions, changes in the surface abundances prior to 
the asymptotic giant branch (AGB) stage are only ex- 
pected to be due to convective dilution during the first 
dredge-up phase. In low mass red giant branch (RGB) 
stars, the convective envelope reaches only regions where 
^^C was processed to ^'^C and ^''N. Consequently, the car- 
bon isotopic ratio declines (from 90 to about 20-30), the 
carbon abundance drops (by about 30%) and nitrogen in- 
creases (by about 80%), but oxygen and all heavier el- 
ement abundances remain unchanged. According to the 
standard scenario, the surface abundances then stay unal- 
tered as the convective envelope slowly withdraws outward 
in mass during the end of the RGB evolution. 

However, observational data reveal a different reality. 
The first discrepancy came among the first stellar ^^C/^'^C 
ratios published (Day et al. 1973), when Arcturus was 
found to have ^^C/^'^C = 7.2 ± 1.5. The discrepancy was 
aggravated by the analysis by Lambert & Ries (1981) of 
C, N, and O abundances in a sample of red giants, in- 
cluding Arcturus. In fact, in most of the metal-deficient 
field and globular cluster evolved stars, the observed con- 
version of ^^C to ^'^C and ^''N greatly exceeds the lev- 
els expected from standard stellar models; the ^^C/^'^C 
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ratio even reaches the near-equilibrium value in many 
Pop II RGB stars (Sneden, Pilachowski & Vandenberg 
1986; Smith & Suntzeff 1989; Brown & Wallerstein 1989, 
1992; Gihoy & Brown 1991, henceforth GB91; Brown, 
ct al.l990, henceforth BWO; Bell et al. 1990; Suntzeff & 
Smith 1991; Shetrone et al. 1993, henceforth SSP93; Bri- 
ley et al. 1994, 1997). This problem actually also occurs, 
but to a somewhat lower extent, in evolved stars belong- 
ing to open clusters with turnoff masses lower than 2Mq 
(Gihoy 1989; GB91). 

In addition. Population II evolved stars present other 
chemical anomalies. In halo giants, the lithium abun- 
dance continues to decrease after the completion of the 
first dredge- up (Pilachowski et al. 1993). A continuous 
decline in carbon abundance with increasing stellar lu- 
minosity along the RGB is observed in globular clusters 
such as M92 (Carbon 1982; Langer et al.l986), M3 and 
M13 (Suntzeff 1981), M15 (Trefzger et al.l983), NGC 6397 
(Beh et al.l979, Briley et al.l990), NGC 6752 and M4 
(Suntzeff & Smith 1991). In some globular clusters (M92, 
Pilachowski 1988; M15, Sneden et al.l991; M13, Brown et 
al.l991, Kraft et al.l992; w Ccn, Paltoglou & Norris 1989), 
giants exhibit evidence in their atmospheres for O— >N pro- 
cessed material. In addition to the O versus N anticorrela- 
tion, the existence of Na and Al versus N correlations and 
Na and Al versus O anticorrelations in a large number 
of globular cluster red giants has been clearly confirmed 
(Drake et al.l992, Kraft et al.l992, 1993, Norris & Da 
Costa 1995, Kraft 1994, Kraft et al.l997, Shetrone 1996, 
Zuckeret al.l996). 

These observations suggest that low mass stars un- 
dergo a non-standard mixing process which adds to the 
first dredge-up and modifies the surface abundances. The 
first and only evidence on the evolutionary state at which 
this non-standard mixing actually becomes effective came 
from observations of the C/N and ^^C/^^C ratios (respec- 
tively by Brown 1987 and by GB91) in evolved stars of 
M67. In this old open cluster, which has a turnoff mass 
of 1.2 Mq (VandenBerg 1985, Meynet et al.l993) and a 
value of [Fe/II]=0.0±0.1, the observational and theoretical 
first dredge-up appear to be in complete agreement (Char- 
bonnel 1994, henceforth C94). Indeed, the i^c/i^C ratio 
changes at the base of the giant branch from a value >40 to 
21-25 between Mv ~ -h3.5 and Mv « +2.5; C/N changes 
by about a factor of 4 in this same interval. The onset 
of the dredge-up occurs at a luminosity in accord with 
standard theory, and the observed post-dilution ratios are 
in very good agreement with the standard predictions at 
this point. However, in M67 an additional mixing event 
begins at about Mv w -1-0.5 where ^^C/^^C drops from 21- 
25 to 11-15. This luminosity actually corresponds to the 
so-called "bump" in the RGB luminosity function (Fusi- 
Pecci et al.l990), i.e., to the evolutionary point where the 
hydrogen burning shell crosses the chemical discontinu- 
ity created by the convective envelope at its maximum 
extent. As pointed out by C94, this observational fact 



strongly suggests that prior to this evolutionary point the 
mean molecular weight gradient created during the first 
dredge-up acts as a barrier to any mixing below the con- 
vective envelope. After this point, however, the gradient 
of molecular weight above the hydrogen-burning shell is 
much lower, and extra-mixing is free to act. 

Since the observations in M67, there have been no new 
sets of abundance data in any cluster, either globular or 
galactic, to establish the evolutionary point at which the 
extra-mixing begins to occur. With 4-metcr class tele- 
scopes it is almost impossible to observe red giants in 
globular clusters at the luminosities at which ^^C/^^C is 
expected to drop from near 20 to less than 10 because 
the molecular bands become weaker with increasing tem- 
perature and gravity. Furthermore, the attainable signal- 
to-noise of the spectra diminishes with decreasing stellar 
luminosity. It is crucial, however, for our understanding of 
the nature of the process to compare conditions in stars 
of different populations and metallicitics. 

To this end, we address here the problem of the onset 
of the extra-mixing, and of its possible metallicity depen- 
dance, by placing in an evolutionary sequence five field 
giants with moderate metal deficiencies and two giants 
of the globular cluster 47 Tuc for which we have assem- 
bled accurate abundance data for Li, C and N as well 
as for the ^^C/^'^C ratio. Our observational data are pre- 
sented in §2. We enlarge our sample by including in our 
discussion the old disk giants with ^^C/^^C ratios derived 
by SSP93. Very precise absolute magnitudes were recom- 
puted from HIPPARCOS parallaxes for all the stars we 
consider. As discussed in §3, the sequence allows us to con- 
firm that the discrepancy between observations of mixing- 
sensitive species and standard theory really appears at the 
luminosity-function bump (LFB) on the RGB, but not at 
lower luminosities. Since our sample is moderately metal- 
deficient, we can study the metallicity-dependence of the 
extra-mixing by comparing our results with the observa- 
tions in M67. In §4, we investigate the nature of the fac- 
tor that stabilizes the star against extra-mixing before it 
reaches the luminosity-function bump, and discuss the role 
of molecular- weight (/i) barriers. From simple considera- 
tions, we derive from our data the "observational" critical 
/i-gradient, (Vln/i)c.o6s, that shields the central regions of 
a star from extra-mixing. If we assume that extra-mixing 
is free to act down to the layers defined by the same value 
of (V In /i)c,ofcs whatever the stellar mass and metallicity, 
then ^^C/^'^C ratios can be explained both for Pop I and 
II giants. Finally, we discuss the implications for the na- 
ture of the extra- mixing mechanism itself on the RGB, 
and bring clues that it may be related to rotation-induced 
processes. 

2. Observational data for metal poor stars 

Some of the data that we employ here have already been 
published. We concentrate on stars of modest metal defi- 
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Table 1. Observed Stars 



HD 


Star 


V- K 


T.// 


Mv 


Mi,o; 


logg 


[Fe/H] 


37160 


4? Ori 


2.41 


4720 


+1.33± 0.07 


+0.89 


2.45 


-0.70 


3546 


£ And 


2.13 


4925 


+0.77± 0.09 


+0.43 


2.34 


-0.66 


89484 


7Leo A 


— 


4440 


-0.32±0.07 


-0.9 


1.62 


-0.46 


89485 


7Leo B 


— 


4900 


+0.87±0.07 


+0.5 


2.36 


-0.42 


124897 


Arcturus 


3.00 


4330 


-0.30±0.02 


-0.98 


1.55 


-0.62 


— 


47 Tuc 3501 


3.45 


3975 


-1.1 


-1.9 


0.9 


-0.65 


— 


47 Tuc 4418 


3.50 


3975 


-1.1 


-1.9 


0.9 


-0.65 



° The stars were selected so as to be as close to the evolutionary tracks delineated by the color-magnitude diagram of 47 Tuc 
as their parallaxes (ground based when we selected the stars) and [Fe/H] values from the literature permitted 



Table 2. Comparison of our abundances (upper analysis) with those of Cottrell and Sneden (1986) and Shetrone et al.(1993) 
(lower analysis) of e And and <^^Ori 



Star 


Tcff 


logg 


[Fe/H] 


[C/Fc] 


[N/Fe] 


[0/Fe] 


12^/13^ 


loge(Li) 


£ And 


4925 


2.55 


-0.45 


-0.29 


+0.24 


-0.06 


7 


< -0.2 




4750"'* 


2.0"'* 


_Qr,a,6 


-0.15* 


-0.10* 


+0.2* 


7.5" 


0.0* 


<^2 Ori 


4720 


2.85 


-0.43 


-0.03 


+0.08 


+0.15 


>20 


+0.3 




4600 


2.4 


-0.7 


0.00 


-0.20 


+0.15 


>40 





" From Shetrone et al. (1993) 

* From Cottrell and Sneden (1986) 



ciency. Five bright metal-deficient field stars are 0^ Ori, e 
And, 7 Leo B, 7 Leo A and Arcturus. Their Mboi values 
range from near +0.9 to —1.0 so they form an evolutionary 
sequence that starts from fainter absolute magnitude than 
can be reached with high spectral resolution in globular 
clusters. As a continuation of our "evolutionary sequence" 
we include two stars from the globular cluster 47 Tucanae 
whose properties indicate that they are first giant branch 
(i.e. not AGB) stars. The basic properties of our sample 
stars are shown in Table 1. 

For the field stars we derive the gravity values as- 
suming a mass of 0.8 Mq, the HIPPARCOS parallaxes. 
Toft's derived from the V — K color-Teff relation (Ridgway 
et al.l980, DiBenedetto 1993, Dyck et al.l996), and the 
bolomctric corrections of Buser and Kurucz (1992). This 
method is probably superior to spectrocopic determina- 
tion of gravities using ionization balance; a test using the 
HIPPARCOS parallaxes indicates the spectroscopic grav- 
ities often are too low by a factor of two (Nissen & H0g 
1997). The visual binary stars in 7 Leo do not have sep- 
arate colors; instead, we adjust the spectroscopic Teff's of 
Lambert & Ries (1981) downward by 210K, the mean off- 
set between their temperature scale and the V — K scale. 
We note that the 47 Tuc stars are probably not on the 
same absolute magnitude scale as the HIPPARCOS re- 
sults, since they depend upon a cluster distance modulus 
which probably needs revision in light of the HIPPARCOS 



results for the subdwarfs; see, e.g., Reid (1997). For our 
purposes, the precise value of the absolute magnitudes of 
the two globular cluster stars is not important: it matters 
here only that they are ~ 1 mag more luminous than our 
most luminous field stars, tl 

The abundances we use have been rederived with a 
number of improvements. For Arcturus and the 7 Leo pair 
we have combined the equivalent widths of Lambert & 
Ries (1981) and Ries (1981) with the empirical photomet- 
ric Tcff scale and the most recent and almost certainly 
most accurate value of the CN dissociation potential D° 
= 7.65 eV (Bauschficher et al.l988). For 47 Tuc we have 
used the data of BWO with two changes. The analysis 
by Hesser et al.(1987) of the color-magnitude diagram of 



^ Both Arcturus and Gamma Leo A lie above the horizontal 
branch and hence could be post-He-flash stars. However their 
colors place them on the first ascent giant branch. A compar- 
ison of the colors we have quoted shows that Arcturus falls 
exactly on the first ascent giant branch of 47 Tuc (Hesser et 
al. 1987). For Gamma Leo the measured B-V color pertains 
to both stars and must be deconvolved to obtain the color of 
star A. To derive the observed B-V color of 1.15 mag, we com- 
bined the color of Arcturus, B-V=1.23 with that of the G8III 
secondary, 0.95 mag, using the observed 1.2 visual magnitude 
difference, to derive the observed composite color of 1.15 mag. 
This confirms that the color of Gamma Leo A is about the 
same as that of Arcturus and hence that Gamma Leo A on the 
first ascent giant branch. 
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Table 3. Abundance Results 



Star 


12^ 


13c 


N 


O 


(C + N)" 


(C + N + O)*^ 


12^/13^ 


log 6(Li) 


(j>^ Ori 


7.9 


<6.6 


7.8 


8.55 


8.17 


8.7 


>20'= 


0.3 


e And 


7.65 


6.8 


7.95 


8.3 


8.15 


8.53 


7d 


< -0.20" 


7 Leo B 


7.77 


6.8 


7.87 


8.45 


8.14 


8.62 


9.5 


< +0.3-'^ 


7 Leo A 


7.8 


7.0 


7.9 


8.45 


8.18 


8.64 


6.5 


< 0.0-'^ 


Arcturus 


7.9 


7.0 


7.7 


8.55 


8.15 


8.7 


7 


< -1.5» 


47 Tuc 3501 


7.7 


6.85 


8.25 


8.45 


8.37 


8.7 


8 


h 


47 Tuc 4418 


7.6 


6.75 


8.25 


8.5 


8.35 


8.7 


7 


h 


0' 


8.6 


6.6 


8.0 


8.9 


8.7 


9.1 


90 


1.15 



This is ^^C + i^C + N 

This is i^C + ^^C + N + O 

Shetrone et al.l993 find ^^C/^-'^C > 40 

Shetrone et al.l993 find ^^C/^'^'C = 7.5 

Shetrone et al.l993 find Li = 0.0 

Heifer & Wallerstein 1968 

Lambert et al.l980 

not present, no upper limit computed 

photospheric value, Anders & Grevesse 1989, Grevesse et al.l991 



47 Tuc indicate an iron abundance for 47 Tuc higher by 
about 0.2 dex as compared with the BWO value. Hence, 
we have accepted [Fc/H] = —0.65 rather than —0.85 for 
the 47 Tuc stars and have recalculated the CNO abun- 
dances using models with the new [Fe/H] values. For the 
12(-;^13q j.g^^jQ jj^ 47 rp^(, g^-jjj. 4428, we have adopted the 

second solution in Table 2 of Bell et al.(199G; multiple so- 
lutions for ^^C/^^C and carbon abundance are given for 
different values of the stellar microturbulence parameter) 
because the corresponding derived carbon abundance is 
closer to that of BWO (in which the microturbulence is 
derived as part of the analysis of the high resolution spec- 
tra in the CH region). 

To extend our evolutionary sequence to lower luminos- 
ity on the giant branch we have obtained new high reso- 
lution spectra of two field stars (jP Ori and e And, whose 
metallicities are similar to those of the other stars. The 
new spectra were obtained at the Dominion Astrophysical 
Observatory with the 1.2-meter telescope, coude spectro- 
graph, long camera and Reticon detector. This combina- 
tion yields a resolving power of 35,000, and the exposures 
were timed for a signal-to-noise of about 200. Our spectra 
cover the C2 Swan band at 5635 A, the [O I] line at 6300 
A, the Li line at 6700 A and the 2-0 red system CN band 
at 8000 A, thereby providing abundances of Li, C, N, O as 
well as iron and the ratio ^^C/^'^C. These stars were an- 
alyzed using the same models and gf-values as were used 
in the other analyses except that the C2 Swan bands re- 
placed the violet CH bands as the source of the carbon 
abundance. To analyze the C2 Swan band at 5635 A it 
was necessary to synthesize the spectrum over an inter- 
val from 5626-5640 A. The spectroscopic data for the C2 



Swan bands is that reviewed by Grevesse et al.(1991) in 
their re-evaluation of the solar C abundance. 

Both £ And and cfP' Ori have been analyzed by others, 
most recently with modern model atmospheres and digital 
spectra by Cottrell & Sneden (1986) and subsequently by 
SSP93. These two stars provide an overlap with their large 
sample of stars of generally lower luminosity. We compare 
their data with ours in Table 2. 

Most of the differences between our abundances and 
those of the Texas group can be understood as due to 
differences in effective temperatures and gravities. Their 
lower effective temperatures are responsible for their lower 
iron abundances and nitrogen abundances because of 
lower opacities and the dissociation equilibrium of the CN 
molecule. The low masses derived by SSP93 (Table 6) and 
of Cottrell & Sneden (1986, Table 6) demonstrate that 
their gravities are too low by roughly a factor 3 which is 
the difference between their and our gravity values for the 
same stars. 

Our final adopted abundances are shown in Table 3. 

It is important to understand the uncertainties in the 
individual abundances in Table 3. The iron abundances 
are almost independent of the CNO abundances but not 
entirely so. Since many lines of Fe I and several lines of 
Fe II are easily observed the measured line strengths are 
not an important source of error. The uncertainty is surely 
dominated by the uncertainty in the atmospheric model, 
especially at small optical depth where mechanical en- 
ergy input and backwarming, such as electron conduction 
from the chromosphere, can effect the temperature near 
the boundary, especially the temperature minimum (Ku- 
rucz 1996). Hence, we are dealing with an uncertain func- 
tion, T(t), not just an uncertain parameter. This problem 
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was evaluated carefully by Leep et al.(1987) in connection 
with the metallicity of M71 and the derivation of solar f- 
values using either the Holweger-Miiller (1974) model or 
the Bell et al.(1976) model. The conclusion of that pa- 
per was that no [Fe/H] value for a red giant can be more 
accurate than ±0.2 and we accept that here. The influ- 
ence of CNO abundances on the [Fe/H] value comes in 
through the opacity of CO and CN which effects the at- 
mospheric models. The infrared CO bands are important 
for the cooling and structure of the outermost layers of an 
atmosphere of a late-type star, but this effect is less im- 
portant for stars as warm as those in this study. The CN 
lines can contribute substantially to the atmosphere struc- 
ture of red giants because the CN red system blankets a 
rather large wavelength range near the flux peak, but CN 
diminishes in importance quickly as the metallicity of a 
star is decreased because the CN column density falls ap- 
proximately as metallicity squared while the continuous 
opacity falls approximately as metallicity. 

The uncertainty in the oxygen abundance is mostly 
determined by the actual measurement of the equivalent 
width of the only useful oxygen hne, A6300, in the globular 
cluster stars, and the second line at A6363 which is barely 
detectable in the bright field stars. Since the oxygen line 
is formed over a wide range of optical depth the derived 
abundance is not very sensitive to the boundary tempera- 
ture. The low carbon abundances in these stars means that 
the uncertainty in the C abundance has little effect on the 
solution for the O abundance despite the formation of CO 
in the atmosphere. For oxygen the uncertainties in 0/Fe 
are surely ±0.1 dcx or perhaps a little lower for Arcturus, 
for which the measurements are by far the most accurate. 
Hence, the somewhat higher 0/Fe value for Arcturus, as 
compared with the other stars in Table 3, is probably real. 

The formation of CO depletes carbon and an error of 
0.10 dex in oxygen will cause an error of 0.05 dex in the 
carbon abundance. Fortunately the CH and C2 lines are 
not formed at a small optical depth (as they are in the sun) 
but rather below the levels at which carbon is depleted by 
CO formation, so the uncertainty in the boundary temper- 
ature is not important. Considering both the uncertainty 
in the oxygen abundance, the actual measurements of CH, 
and the presence of some saturation in the CH lines, an 
uncertainty of 0.15 dex is possible for the ^^C abundances. 

The ^^C/^'^C ratio is determined largely by the clump 
of ^■'CN lines near 8005 A. This was illustrated in Figure 
3 of BWO: for that star, M13 1-48, a ratio of 6 was derived 
and the figure shows that a range from 5 to 8 is possible. 
Scaling that evaluation to Table 2 of this paper indicates 
that an observed ratio of 7 for three stars could lie between 
about 5.5 and 9; the pure equilibrium ratio of 3.5-4.0 is 
pretty well excluded as are values of 10 or above. 

The nitrogen abundance, derived from CN lines, is 
less dependent on line measurement uncertainties because 
many (^15-20) lines are available in our spectra. However, 
it does depend on the abundances of C and O through the 



dissociation equilibrium of CO, CH, and CN, as well as 
the formation of N2. Hence the uncertainties in the C and 
O abundances come into play as well as the uncertainty 
in the CN dissociation potential. For CN the difference 
between the theoretical value, 7.65 eV (Bauschlicher et 
al.l988) and the recent experimental value of 7.77 ± .05 
eV (Costes et al.l990) is disturbing. We have used the 
7.65 eV value. If the 7.77 eV value is correct our N abun- 
dances must be lowered by 0.2 dex, which lowers our C±N 
values by a somewhat smaller amount. The reduction in 
C±N±0 introduced by such a reduction in N is .05 dex 
for Arcturus and the 7 Leo pair, and .08 dex for the 47 
Tuc stars. 

An error in either the dissociation potential of CN or 
the /-value of the CN red system 2-0 band would effect all 
of our nitrogen abundances equally. Both of these quan- 
tities are not as well determined as might be hoped: the 
dissociation energy of CN seems to be uncertain by about 
0.2 eV, and different sources for the f- values disagree by as 
much as 30(see the discussion in Bauschlicher et al.l988). 
The values used here for these quantities, 7.65 eV for Dq 
and 8.4 X 10^^ for /20, are those used in (for Dq) or derived 
from (/20) the analysis of CN red system lines in the solar 
spectrum by Sneden & Lambert (1982), so by construction 
they yield the solar nitrogen abundance. A second source 
of uncertainty in the N abundance is introduced by the 
uncertainty in the C abundance. For measured CN lines 
any error in the C abundance introduces an error of com- 
parable size and opposite sign in the N abundance. Since 
an error of as much as 0.15 dex in C is possible, an error 
of the same size with reversed sign is possible for N. 

The ^^C/"'^^C ratio also has some observational uncer- 
tainty, almost exclusively due to the difference in strength 
between the ^2CN and ^3CN lines, so that the microturbu- 
lent velocity enters into the isotope ratio solution. This un- 
certainty is aggravated as the ^ 2CN lines get stronger and 
as the isotope ratio increases. For the stars in this sample, 
the -'^2CN lines are quite weak: for 7 Leo A, which has the 
strongest CN lines in the sample, the reduced equivalent 
widths are in the range —5.5 < log{W/X) < —5.1, indi- 
cating that saturation (and hence microturbulence-related 
errors in the isotope ratio) are minimal. Our results are 
summarized in Figures 1 and 3, and discussed below. 

3. The onset of extra-mixing on the observed evo- 
lutionary sequence 

Our observational results are summarized by referring to 
Figures 1 to 3. We show in Figure 1 the observed depen- 
dence of ^^C/^^C ratio on Mboi- This sequence for our 
sample of stars with moderate metal deficiencies reveals 
a behavior of the surface ^^C/^^C ratio very similar to 
the one previously observed at solar metallicity in M67. 
Three main features appear: (i) The least luminous star, 
(fp Ori, shows a depletion of Li by 2.7 dex below the level 
accepted for young stars and a limit on the ^^C/^'^C ratio 
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Fig. 1. The ratio of ^^C/^'^C is plotted against the bolometric 
magnitude of the program stars. Uncertainties of ±2 for the 
isotope ratios are estimated for each star except for Arcturus, 
for which ±1 is estimated. The absolute magnitudes are derived 
from the HIPPARCOS parallaxes. Note that 47 Tuc 3501 and 
4418 have both the same luminosity and carbon isotopic ratio. 
The horizontal line indicates the position of the bump in the 
luminosity function of 47 Tuc (King et al.l985). 



Fig. 2. Carbon isotopic ratios for the old disk giants by Cottrell 
& Sneden (1986; crosses) and Shetrone et al.(1993; open circles 
and triangles)). Our program stars are also shown (black circles 
and triangle). Observations for the same stars are connected 
by broken lines. The absolute magnitudes are derived from 
the HIPPARCOS parallaxes. The horizontal line indicates the 
position of the bump in the luminosity function of 47 Tuc (King 
et al.l985). 



in agreement with standard predictions for dilution El This 
indicates that the standard theoretical main sequence pro- 
files of ^^C and ^^C match the real chemical profiles, and 
that the extra-mixing is only efficient after the comple- 
tion of the first dredge- up (see C94). (ii) Then between 
Mboi = -1-0.9 and -1-0.5, the observed isotopic ratio drops 
to values near 7, well below the standard predicted post- 
dilution ratio. It is interesting to note that in the RGB 
luminosity function of 47 Tuc, King, et al.(1985) local- 
ize the bump at Mboi = -1-1.05 ± 0.2, i.e. precisely in the 
region where the disagreement between standard predic- 
tions and observations of the carbon isotopic ratio appears 
in our sample. This confirms that the extra-mixing which 
leads to very low ^^C/^^C ratios in low-mass and metal- 
deficient evolved stars becomes efficient exactly when the 
hydrogen-burning shell crosses the chemical discontinuity 
created by the outward moving convective envelope, (iii) 



^ (All the standard models computed by different authors 
give very similar result for the post-dredge up value of 
12C/13C; see for example the compilation by Wasserburg et 
al. 1995. Boothroyd & Sackman (1997) showed that different 
initial carbon isotopic ratio do not affect significantly the post- 
dredge up result, whatever the metallicity of the models.) 



Finally, from Mboi = +0.4 to —2, there is no further change 
in the ^^C/i^C ratio. 

The "abrupt" change in ^^C/^^C ratio shortly after 
the luminosity-function bump occurs both at solar and 
at lower metallicities, although the final -'^^C/^'^C ratio is 
lower in the more metal-poor stars. It is worth noting that 
the change in ^^C/^'^C ratio is "abrupt" only in terms of 
stellar luminosity. The LFB is caused by a slower rate 
of evolution for stars at this luminosity (that is, dL/dt 
is smaller there) owing to the H-burning shell contacting 
the H-rich, previously mixed, zone; of the time stars take 
to evolve from the beginning of the LFB to the tip of 
the giant branch, 15 — 20% of it is spent in the LFB. 
Consequently, the change in the ^^C/^^C ratio may not 
be so abrupt in terms of time. 

In order to enlarge our sample, we recomputed the 
absolute magnitudes from the HIPPARCOS parallaxes 
for the old disk giants with ^^C/^'^C ratios derived by 
SSP93. These stars also have moderate metal deficiencies 
(-1.0<[Fc/H]<-0.3) . Due to the good luminosity determi- 
nation, the complete data, which are presented in Figure 
2, can be reasonably viewed as an evolutionary sequence. 
We also show the ^^C/^'^C ratios previously derived by 
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CottrcU & Sneden (1986) for the stars which belong to 
the SSP93 sample. 

The SSP93 isotope ratios tend to be systematically 
higher than the Cottrell & Sneden values for stars with 
high ^"^C/^^C, but this is easily understood as being due 
to SSP93's superior spectral resolution. The higher spec- 
tral resolution allows detection and better measurement 
of the very weak ^'^CN features. This highlights a gen- 
eral problem in the determination of ^^C/^^C ratios: three 
effects operate to make ^^CN detection difficult for low- 
luminosity stars. First, at lower luminosities the ^^CN pro- 
portion is lower due to their unmixed state. Second, model 
atmosphere effects operate so that higher-gravity, lower- 
luminosity stars have weaker CN lines even for constant 
chemical composition in the atmosphere. Third, the lower 
stellar luminosity makes it more difficult to obtain the 
requisite S/N to detect the weak lines. 
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Fig. 3. Log '^^C/'^^N against Mboi for the program stars. The 
uncertainty in log ^'^C/^^N is estimated to be 0.2 dex except for 
Arcturus, for which it is probably 0.1 dex. The horizontal line 
indicates the position of the bump in the luminosity function 
of 47Tuc (Kinget al.l985). 



Despite the uncertainties, it is clear that in the present 
complete sample, no star with Mboi > +1-0 shows a car- 
bon isotopic ratio lower than the one predicted by stan- 
dard dilution models. Figure 2 is strong evidence that ex- 
tra mixing begins at a luminosity of about Mboi = +1- 
This limit in luminosity corresponds to the approximative 
value for Mboi (bump) observed in clusters of the same 
[Fe/H] as the present stars (Fusi Pecci et al.l990). The 



cutoff region appears broader here than in our restricted 
sample, mainly due to a probable scatter in the stellar 
masses in SSP93 data. Indeed, a O.IM© difference corre- 
sponds approximatively to a shift of 0.2 in Mboi for the 
position of the bump (Fusi Pecci et al.l990). 

We show the observed ^^C/^'*N ratios plotted against 
Mboi in Figure 3. (f^ Ori presents a ratio consistent with 
standard post-dredge-up predictions. The ^^C/'^^'N ratio 
decreases simultaneously with the ^^C/^^C ratio between 
Mboi = -1-1.0 and -1-0.5 by about 0.3 dex, reveahng that the 
extra-mixing process also slightly affects the N abundance. 
Arcturus, for which the data is by far the best due to its 
brightness, stands however above the curve for the other 
stars (which really means above 7 Leo A, at the same 
luminosity, by 0.3 dex). u This difference was recognized 
easily 15 years ago by Lambert & Ries (1981), who found 
a difference of 0.47 dex. The 47 Tuc stars, at Mboi = -1-9, 
present ^^C/^'*N ratios lower than those in the less evolved 
objects. 

In the M67 data (Brown 1987), a change in the C/N 
ratio at and above the luminosity of the RGB bump is 
not clear. This is probably due to errors, both random 
and systematic, in the abundances derived from the syn- 
thetic spectrum fits, and aggravated by the sparseness of 
the cluster RGB. A systematic offset in the nitrogen abun- 
dance was found in that study, and an attempt was made 
to calibrate this out in a post hoc way. The result, how- 
ever, is that any slow trend in C/N with luminosity in the 
Brown (1987) data must be interpreted with caution. It is 
more secure, however, to note that the total change in the 
C/N ratio between stars at the base of the giant branch 
and those at the RGB tip and in the post-helium core flash 
clump is larger than predicted in standard models, so that 
the most luminous M67 giants again have also engaged in 
some extra mixing. 

In our sample, the Li abundance drops between <jp 
Ori and e And (in terms of luminosity) after which it 
is below detectability. It is difficult however to argue for 
the luminosity onset of extra-mixing from this data. In- 
deed, lithium destruction on the main sequence may lead 
to a dispersion of the abundances on the red giant branch 
which could explain the difference between 0^ Ori and e 
And. In any case, the lithium abundances in the more 
evolved stars from our sample clearly indicates that an 
extra-mixing mechanism transports lithium from the con- 
vective envelope down to the region where it is destroyed 
by proton capture after the end of the dilution phase. In a 
large sample of evolved halo stars, Pilachowski, Sneden & 
Booth (1993) also showed that the lithium abundance con- 
tinues to decrease after the completion of the first dredge- 
up. 



'^ Aoki and Tsuji (1997) have reanalysed the N abundance 
in Arcturus (and many cooler stars) and find log(^'^C/^*N) = 
0.09, which is within our uncertainty and somewhat closer to 
the 12C/14N ratio in Gamma Leo A. 
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4. Determination of a critical /x-gradient from the 
observed carbon isotopic ratios 

Recently, different groups have simulated extra-mixing be- 
tween the base of the convective envelope and the hy- 
drogen burning shell in order to reproduce the CNO 
abundances in RGB stars. Denissenkov & Weiss (1995) 
modeled this deep mixing by adjusting both the mixing 
depth and rate in their diffusion procedure. Wasserburg et 
al.(1995) and Boothroyd & Sackman (1997) used an ad- 
hoc "conveyor-belt" circulation model, where the depth 
of the extra-mixed region is related to a parametrized 
temperature difference up to the bottom of the hydrogen- 
burning shell. On the other hand, other authors attempted 
to relate the extra-mixing with physical processes, among 
which rotation seems to be the most promising. Sweigart 
& Mengel (1979) suggested that meridional circulation 
on the RGB could lead to the low ^^C/^^C observed in 
field giants. Taking into account the interaction between 
meridional circulation and turbulence induced by rotation 
in stars, Charbonnel (1995) showed that rotation-induced 
mixing can account for the observed behavior of carbon 
isotopic ratios and for the Li abundances in Population II 
low mass giants. 

In these different approaches, the common underlying 
question is the determination of the extension of the re- 
gion where extra-mixing occurs, and more precisely the 
nature of the stabilizing factor. An investigation into the 
mechanism of extra mixing must include an understand- 
ing of why the mechanism does not operate at luminosities 
below that of the RGB bump, and what is the nature of 
the factor that stabilizes the star against extra mixing be- 
fore it reaches this point. When suggesting that the extra- 
mixing on the RGB could be related to rotation, Sweigart 
& Mengel (1979) and Charbonnel (1994, 1995) underhned 
the inhibiting effect of molecular weight (or /i) barriers. 
Indeed, as argued by Mestel (1957), gradients of molecu- 
lar weight tend to restrain the circulation and to stabilize 
the mixing processes related to rotation. We will show now 
how the observations gathered in this paper allow us to 
constrain, from simple considerations, the conditions in 
which a /i-barrier shields the central regions of a star from 
extra-mixing. 

We have computed some standard stellar models for 
exploration of these questions. The evolutionary sequences 
were followed from the pre-main sequence up to the RGB 
tip with the Toulouse- Geneva code (see Charbonnel et 
al.l992). The input microphysics is the same as used in 
Charbonnel et al. (1997). The relative ratios for the heavy 
elements correspond to the mixture by Grevesse & Noels 
(1993). We take the same isotopic ratios as Maeder (1983). 
Neither element segregation or rotation-induced mixing 
are included in the present computations. 

The gradients of molecular weight ( V In y^ = d In n/dr) 
in a standard stellar model of 0.9Mo with [Fe/H] = —0.45, 
typical of the present observed sample, are presented in 



0.9Mq, [Fe/H]--0.45 
end of main sequence 
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Fig. 4. Molecular weight gradient (in log scale, bold dashed 
line, left axis) and composition profiles (right axis) in a O.QM©, 
[Fe/H] =-0.45 stellar model, at the end of the main sequence. 
The mass fractions are multiplied by 100 for "^He and ^''N, by 
1000 for "C; the ratio ^^C/^^C is divided by 100. The thick 
cross corresponds to (Vto/it)c, the critical molecular weight 
gradient obtained with the prescription by Huppert & Spiegel 
(1977) discussed in §5 



Figures 4, 5 and 6 at three different evolutionary stages: 
at the end of the main sequence, after the completion of 
the first dredge-up, and shortly after the RGB bump. We 
also show the abundance profiles for H, "^He, ^^C and ^''N, 
and the ^^C/i^C ratio. 

While on the main sequence, nuclear burning leads to 
an increase of the mean molecular weight in the central re- 
gions of the star (see Figure 4) . During the first dredge-up 
phase, the convective envelope homogenizes the star down 
to very deep regions, and builds a very steep gradient of 
molecular weight at the point of its maximum penetration 
(Vln^ ~ 10"^ at Mr/M ~ 0.282 in Figure 5). Further on, 
as the star keeps ascending the RGB, the hydrogen burn- 
ing shell becomes thinner and moves outwards in the mass 
scale, while the convective envelope retreats. As discussed 
in §3, our observations provide strong evidence that the 
extra-mixing is inhibited until the star reaches the LFB. 
The shape of V In /i just after this evolutionary point, i.e. 
after the encounter of the hydrogen burning shell with the 
previously mixed region, is shown in Figure 6. Below the 
base of the convective envelope V In /i is much lower at 
this time. 

Our observations provide a precise clue on the exten- 
sion of the extra- mixed region down to the nuclear burning 



Charbonnel, Wallerstein and Brown 



-10 



O 



-15 



-20 



end of 1st drcdgc-up 



/^ 



0.8 



\ - 0.6 




0.4 



0.2 



0.24 



0.25 



0.26 0.27 

M /M 



0.28 



0.29 



Fig. 5. Same as Fig. 4, after the completion of the first 
dredge-up, in the H-burning shell region. The base of the con- 
vection zone is located at Mr/M~0.283. 
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Fig. 6. Same as Fig. 4, after the RGB bump, in the H-burning 
shell region. The base of the convection zone is located at 
Mr/M~0.313. 



layers. The low ^^C/^^C ratios we see in our sample (~ 7) 
can indeed be reached if the extra-mixing develops down 
to Mr/M~0.294 in our 0.9Mq model (see Figure 6). At 
this point, Vln/i ~ 1.5 x 10~^^. We now assign this value 
to (V ln/i)c,o6s, that we define as the "observational" crit- 
ical /u-gradient which appears to shield the central regions 
of the star from extra-mixing. In our 0.9Mq model, the 
temperature difference A log T between the bottom of the 
hydrogen burning shell and this point is equal to 0.26. This 
value for A log T is in agreement with the one Boothroyd 
& Sackman (1997) have to impose in their parametrized 
mixing model in order to match the observed ^^C/^'^C ra- 
tios at solar metallicity. 

Let us check now to what region (V In /i)c,o6s — 1-5 x 
10~^^ corresponds, in terms of expected ^^C/^^C, at differ- 
ent evolutionary stages, and for stars of different metallic- 
ities. We see in Figure 4 that while the star is on the main 
sequence, the real V In /x becomes higher than (V In fJ,)c,obs 
in the very external part of the star, around M^/M equal 
to 0.75. If we assume that such a /x-barrier can not be 
penetrated, then on the main sequence no extra-mixing is 
expected to occur in the stellar region of energy produc- 
tion. This explains the observations of the carbon isotopic 
ratio in (f>^ Ori, in the less luminous stars of SSP93's sam- 
ple, and in the less evolved RGB stars of M67 (GB91), 
where we see a perfect agreement with standard theoret- 
ical dilution. This is also in agreement with hclioseismo- 
logical constraints. Indeed, solar structure computations 
bring important information on the critical /i-gradients 
that limit the mild mixing which is necessary to explain 
the lithium depletion in the Sun. Richard et al.(1996) 
and Richard & Vauclair (1997) showed that the best solar 
models (for what concerns helioseismological comparison 
and agreement with Li and Be observations) are those 
including both element segregation and rotation-induced 
mixing, where the mixing is cut-off when the /x-gradient 
becomes > to 1.5 - 4 xlO^^^. It is highly satisfactory to 
get the same value for (V In /i)c,ofcs from two completely 
different observational constraints. 

The gradients of molecular weight shortly after the 
RGB bump in a model of O.9M0 with [Fe/H] = -2 (typi- 
cal of a globular cluster giant) and in a model of 1.25M0 
at solar metallicity (typical of a M67 giant star) are pre- 
sented respectively in Figures 7 and 8. If we assume that 
the extra-mixing reaches down to the same (V In /i)c,ofcs 
whatever the star, then the ^^C/^'^C is expected to de- 
crease down to its equilibrium value for the globular clus- 
ter giant, and down to ~ 12 in the open cluster giant. 
These values, and the dependency with metallicity, are ex- 
actly as observed in the giants in open clusters by Gilroy 
(1989) and are in agreement with the other observations 
of field stars and globular cluster giants. They are inde- 
pendent of any modeling for the extra-mixing. It appears 
then that the value we derived empirically for (V In /i)c,o6s 
is "universal" , in the sense that an unique value is suffi- 
cient to describe different constraints. 
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Fig. 7. Same as Fig. 6 in tlie O.9M0, [Fe/H]=-2 stellar model. 
Tlie base of the convection zone is located at Mr/M~0.387. 
The mass fractions are multiplied by 100 for He, and by 50000 
and 5000 respectively for ^^C and "N; the ratio ^^C/^^C is 
divided by 100 
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Fig. 8. Same as Fig. 6 in the 1.25Mq, [Fe/H]=0 stellar model. 
The base of the convection zone is located at Mr/M~0.218. 
The mass fractions are multiplied by 100, 1000 and 50 respec- 
tively for ^He, ^^C and ^^N; the ratio ^^C/^^C is divided by 
100 



As already shown by C94, if molecular weight barri- 
ers are the inhibiting factor, then no extra-mixing is ex- 
pected to occur in stars which ignite helium in a non- 
degenerate core, i.e. in stars with masses higher than 1.7 
to 2.2M0 (the exact value depends on the metaUicity). In- 
deed, in these more massive stars, the hydrogen burning 
shell does not have the time to reach the regions that have 
been homogeneized during the first dredge-up, so that 
(V In jj) below the convective envelope is always higher 
than (Vln/i)c.o6s- Under these conditions, extra-mixing is 
inhibited in these stars. This is confirmed by the obser- 
vational data in open clusters with turn-off masses higher 
than about 2Mq (Gihoy 1989). 



5. Clues on the nature of the extra-mixing process 
on the rgb 

The results described previously were obtained using ob- 
servational constraints only, and did not rely on any spe- 
cial prescription for the extra-mixing. We can nevertheless 
check whether the value we get for (Vln^)c,o6s is consis- 
tent with what is expected in a very simplified rotation 
framework. Let us follow the suggestion by Huppert & 
Spiegel (1977) according to which meridional currents can 
penetrate into a region with a stable gradient of molec- 

,n(r) 



Nu 



ular weight within a scale height given by h 
where r is the local radius, f2(r) is the angular rotation 
velocity, and N^ the buoyancy frequency due to the /x- 
gradients. We can derive then a critical /i- gradient by 
specifying that h must be a small fraction of r (h = er): 

{S7lnii)c,theor — -^ ^GM( V ' "^ti^re all the quantities are 
computed at the place where the actual V In /i is equal to 
the critical one. 

For each of our 3 models, we computed (V In /i)c,t;ieor 
for a stellar rotation velocity of 2.5 km sec~^, typical 
on the RGB for the stellar masses in consideration (De 
Medeiros et al.l996). The (V In /i)c,theor we obtain in this 
model-dependent approach are shown in Figures 4 to 8 
(thick cross) . They are in perfect agreement with the con- 
straints coming from the observations discussed in §4. This 
result tends to indicate that the extra-mixing process on 
the RGB is related to rotation. 

However, we have to be cautious in our conclusions. 
Indeed, {Vlnn)c,theor is derived here in a simplified frame- 
work. Complete computations with rotation-induced mix- 
ing have to be performed. The transport of chemicals and 
angular momentum should be treated simultaneously to 
properly handle the development of fi-hamcvs (despite ro- 
tational mixing) during stellar evolution. Different stellar 
rotational histories also have to be considered, which may 
result in different efficiencies of the extra-mixing, and lead 
to a dispersion of chemical anomalies from star to star. 
Observational support for this approach may exist from 
the variations of oxygen abundances and rotational ve- 
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locities in blue horizontal branch stars in three globular 
clusters (Peterson et al.l995). 

6. Conclusions 

We have assembled accurate ^^C/^'^C ratios, as well as 
abundance data for C, N, O and Li for five field giants 
with [Fc/H] ~ —0.5 including Arcturus, and for two stars 
in the globular cluster 47 Tuc, and reviewed previously 
published data for a number of other field giants. Using 
HIPPARCOS parallaxes to derive Mboi values, we placed 
these stars into an evolutionary sequence along the as- 
cent of the RGB. The sequence shows that the ^^C/^^C 
ratios drops from ^ 20 to near 7 between Mboi = +0-9 
and +0.5, while Li disappears. The C/N ratio is also af- 
fected, with log (^^C/^^N) diminishing from -1-0.1 to —0.6 
between Mboi = +0.9 and —2. Arcturus stands above the 
curve defined by the other stars by 0.4 dex. These observa- 
tions confirm the existence of an extra- mixing process that 
becomes efficient on the RGB only when low-mass stars 
reach the so-called luminosity-function bump. This result 
is established for the first time for stars with moderate 
metal deficiencies. 

We used the values of the carbon isotopic ratio ob- 
served in our sample to get constraints on the /i-barricrs 
that may shield the central regions of a star from extra- 
mixing. We showed that the same value of the critical 
gradient of molecular weight leads to ^^C/^^C ratios ob- 
served at different metallicities. This observational criti- 
cal /x-gradient is in very good agreement with the critical 
^-gradient which is expected to stabilize meridional cir- 
culation. This result brings strong clues on the nature of 
the extra-mixing which occurs on the RGB, and indicates 
that it is related to rotation. Observational support for 
this approach may exist from the variations of oxygen, 
sodium and aluminum in globular cluster red giants as 
well as the measurement of rotational velcoities in blue 
horizontal branch stars in three globular clusters (Kraft, 
1994; Kraft et al., 1997; Peterson et al. 1995). 
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